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Abstract. Recurrent novae are binaries harboring a very massive white dwarf (WD), 
as massive as the Chandrasekhar mass, because of their short recurrence periods of 
nova outbursts of 10-100 years. Thus, recurrent novae are considered as candidates of 
progenitors of Type la supernovae (SNe la). In fact, the SN la PTF 1 Ikx showed evid- 
ence that its progenitor is a symbiotic recurrent nova. The binary parameters of recur- 
rent novae have been well determined, especially for the ones with frequent outbursts, 
U Sco and RS Oph, which provide useful information on the elementary processes in 
binary evolution toward SNe la. Therefore we use them as testbeds for binary evolution 
models. For example, the original double degenerate (DD) scenario cannot reproduce 
RS Oph type recurrent novae, whereas the new single degenerate (SD) scenario pro- 
posed by Hachisu et al. (1999) naturally can. We review main differences between 
the SD and DD scenarios, especially for their basic processes of binary evolution. We 
also discuss observational support for each physical process. The original DD scen- 
ario is based on the physics in 1980s, whereas the SD scenario on more recent physics 
including the new opacity, mass-growth efficiency of WDs, and optically thick winds 
developed in nova outbursts. 

Keywords : binaries: general - novae, cataclysmic variables - stars; mass-loss, white 
dwarfs - supernovae; general 

1. Introduction 



Nova outbursts are a thermonuclear runaway event on a white dwarf (WD). After the unstable 
hydrogen shell-burning sets in, the envelope greatly expands and the WD becomes very bright 
to reach optical maximum. Structures of such expanded envelopes are subject to the opacity 
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that regulates the energy flux. In the beginning of 1990s opacity table s were revised by the 
OPAL (Iglesias, Rogers & Wilson| [l987; Rogers & Iglesias 1992; Iglesias "&Rogers II 19961) and 
OP dSeaton et al. II 19941) projects. These new opacities showed a prominent peak at a temper- 
ature of logT (K) ~ 5.2 due to iron lines, and brought drastic changes into stellar structures, 
especially when the luminosity is close to the Eddington luminosity. Using the new opacities, 
peop le solved many long-standing problems in pulsation, stellar evolution, and nova outbursts 



(e^. 'Bressan et al II 19931: ISchaller et al]ll992t iKato & Hachisu]ll994t IPrialnik & Kovetzlll995 , 
|G_autschv & Saio II 19961 and references therein). Now a days the OPAL opacities are widely used 
in astronomy. 



Nova theory has been greatly influenced by the new opacities. The strong peak in the opa- 
city accelerates winds during nova outbursts and, as a result, timescales of no va outbursts are 



e.g.: 



Cox & Stewart 



dr astic ally shortened. With t he old opacities (the Los Alamos opacities ^ 

iTOlbl: [Cox & Tabor1 ll976h . we had to assume that novae occur only in very massive WDs 
(MwD ^1.3 Mq) or that the frictional eflrect by a companion star to the WD is very effective in 
ejecting most of the WD envelope in a dynamical timescale. These assumptions were required 
for a relatively short duration of nova outbursts (~ 1 yr). With the new opacities, however, these 
assumptions are not needed anymore. In and after 1990 many numerical calculations were done 
with the OPAL opacities. We have now strong accelerations of winds that blow most of the en- 
velope mass in a relative ly short timescale of ~ 1 yr. As a result, calcul ated nova timescales 
are drastically shortened (iKato & Hachisu 1 1994 ; Prialnik & Kovetz l ll995h . Kato developed an 
optically thick wind theory to follow the extended stage of nova outbursts and her group suc- 
ceeded in r eproducing nova light curves from slow to very fast novae (e.g., Hachisu & Katol 
2006ll2010l) . Now we understood basic properties of novae, such as duration, expansion speed, 
multiwavelength properties, and light curves for all speed classes of novae. 



Type la supemovae (SNe la) are characterized in principle by spectra with strong Si ii but no 
hydrogen lines at maximum light. SNe la play very important roles in astrophysics as a standard 
candle for measuring cosmological distances and as main producers of iron group elements in 
chemical evolution of galaxies. It is commonly agreed that the exploding star is a mass-accreting 
carbon-oxygen (C-hO) W D. For example, the SN la 2011fe indicates its size of the exploding star 



as small as /? < 0.02 Rq (iBloom et al. 11201 2l) . consistent with its WD origin. However, it is not 



clarified yet whether the WD accretes H/He-rich matter from its binary companion (the so-called 
single degener ate (SD) scenario), or two C-i-O WDs merge (the so-caU ed double degenerate (DD) 
scenario) (e.g.. iHillebrandt & Nieme ver II2OOOI i Nomoto et al. II2OO0I for review). It was before 
the advent of the OPAL opacity that Webbink] (Il984 and llben & Tutukovl (Il984 proposed flie 
DD scenario. In this DD scenario, intermediate-mass (3-8 Mq) binaries undergo the first and 
second common envelope phases and finally become a double WD system with a very short 
orbital period of < 3 hr. They can merge within a Hubble time due to orbital energy and angular 
momentum losses by gravitational radiation. If the total mass of a merged object exceeds the 
Chandrasekhar mass, it could explode as an SN la. This is the traditional DD scenario. 



Hachisu. Kato & Nomotd lll996') proposed a new idea, dubbed "accretion wind evolution," 
as a new elementary process in the binary evolution that leads to SN la explosions. This accretion 
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Optically Thick Winds 




Figure 1. Accretion wind evolution in tiie 
SD scenario. Optically thick winds blow from 
a mass-accreting WD when the mass-transfer 
rate from its lobe-filling companion exceeds 
a critical rate, M„. Hydrogen steadily bums 
on the WD at the same rate of (typically 
~ 1 X 10"'' Mq yr' in massive WDs), while 
transferred matter exceeding this burning rate is 
blown in the wind as shown in the figure. 



wind evolution is based on the optically-thick nova wind theory, illustrated in Fig. [T] The WD 
accretes matter from the accretion disk and, at the same time, blows optically thick winds. The 
wind blows as far as the WD accretes matter from the equator at sufficiently high rates. In the 
accretion wind evolution, the strong and fast winds carry mass and angular momentum, so the 
mass transfer can be stabilized. This means that the so-called second common envelope evolution 
does not occur when the ac cretion wind evolution is realized. This changed the traditional binary 
evolution scenarios. Thus. iHachisu et a~ ( 19961) established a new progenitor system of SNe la 
which consists of a mass-accreting WD and a red-giant (RG) donor star In the accretion wind 
evolution, WDs can grow in mass and explode as an SN la. We call this the "accretion-wind 
single degenerat e (awSD) scenario," if i t is necessary to distinguish it from the old SD scenario 
first proposed bv lWhelan & Iben I (Il973h . 



Li & van den Heuvel I (Il997h adopted this idea as a fundamental process of binary evolution 
and found that, in some binaries, the second common envelope evolution does not occur and the 
WD can grow in mass up to the Chandrasekhar limit. These binaries consist mainly of a WD and a 
main-sequence (MS) star at the stage of mass transfer (WD+MS syste ms). This is also a new way 
that lea ds to SN la explosions and dubbed "MS channel." After that, Hachisu. Kato & Nomotd 
1999ah and IHachisu et alJ (Il999bh reformulated basic processes of binary evolution and found 



that the WD-hRG systems (the so-called "symbiotic channel") contribute to the SN la rate as well 
as the WD-hMS systems. 



Recurrent novae are characterized by multiple recorded outbursts of novae. So far, ten objects 
were registered in our Galaxy as listed in Table [T] They should have a very mas sive WD close to 



the Chandrasekhar mass because of the short recurrence time of outbursts (e.g. IHachisu & Kato 



200 Ibt iPrialnik & Kovetz II 19951 see also Fig. |2l). The present evolutionary status of these re- 



current novae can be well understood from the theoretical point of view, if we take into account 
the accretion wind evolution. The binary parameters of several recurrent novae locate just in the 



binary parameter regions of the M S and symbiotic channels of SNe la (e.g., IHachisu & Kato 



|2001b'; Hachisu . Kato.& N omoto 2012b, see also Fig. [3]). Thus we use the binary parameters of 
recurrent novae as testbeds of scenarios on Type la supernova evolution. 
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In the original DD scenario proposed by IWebbink (1984) and Iben & Tutukov (1984), the 



accretion wind evolution was not included, so the binary could evolve to a pair of WDs after 
the second common envelope evolution. Thus, the first essential difference between the SD and 
DD scenarios is in the accretion wind evolution, which was based on the optically thick wind 
theory that was developed after the new opacities appeared in 1990s. In the DD scenario, binaries 
undergo the second common envelope evolution and th e evolutionary p ath is essentially governed 
by the parameters of common envelope evolution (e.g., Webbink The DD scenario is also 

based on other different assumptions from the SD scenario, some of which are based on 1980s 
physics. In the present review, we critically examine these assumptions in light of observations 
of recurrent novae and other related objects. 

Section |2] introduces observational properties of Galactic recurrent novae. Section |3] sum- 
marizes our theoretical understanding of novae, especially some issues closely related to SN la 
progenitors, i.e., stability of accreting WDs and how to determine the WD mass. Section|4]shows 
that these recurrent novae locate at the final stage of evolutionary path to SNe la in the SD scen- 
ario. Section |5] examines several elementary processes that are important in binary evolution. 
Section|6]deals with some objects closely related to the SN la progenitor scenarios. 



2. Recurrent novae: observations 



Recurrent novae are a small subgroup of novae that have multiple recorded outbursts. Their char- 
acteristic properties are summarized as follows. 



1. Recurrence periods between outbursts are as short as ~ 10-100 yrs. 

2. Optical light curves show a rapid decline. 

3. Optical light curves often show a mid-plateau phase. 

4. In the later phase of outbursts, it often becomes a transient supersoft X-ray source (SSS). 
This SSS period almost overlaps with the mid-plateau phase. 

5. No heavy element enhancement is detected in its ejecta. 



Only about ten recurrent novae are known in our Galaxy, but this does not directly suggest a 
small population of recurrent novae. Recurrent nova eruptions develop so fast and their detections 
are not easy especially in their short bright phases. There are potentially a number of recurrent 



novae that are not yet identified as a recurrent nova (Pagno tta et al. 2009). All of the recurrent 



novae are now regarded as a thermonuclear runawa y event owing to instab ility of hydrogen nuc- 
lear burning in a geometrically thin shell on a WD. IWebbink et al. I (119871) categorized recurrent 
novae into three subclasses, in which T CrB and RS Oph were modeled as a pair of an MS and 



an RG, but now it is c lear that they are a binary consisting of a WD and an RG (iSelvelli et al 
19921: 1 Webbink II2OO8I) Also a disk instability model of RS Oph outbursts dXing & Pringle II2009I: 



Alexander et al. I2OI 1|) i s incompatible with many observational indications toward a thermonuc- 



lear runaway event (e.g.. lNelson et al. 
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Table 1. Recurrent novae 



object ^ 


outburst year 


P ^ 

^ rec 


^orb 








(yr) 


(day) 


(day) 


CI Aql 


1917,1941,2000 


24 


0.618 


32 


V394 CrA 


1949,1987 


38 


1.52 


5.2 


V2487 Oph 


1900,1998 


98 


~ 1 


8 


USco 


1863,1906,1936,1945 


8 


1.23 


2.6 




1979,1987,1999,2010 








TCrB 


1866,1946 


80 


228 


6 


RS Oph 


1898,1933,1958,1967 


9 


453.6 


14 




1985,2006 








V745 Sco 


1937,1989 


52 


510 


9 


V3890 Sgr 


1962,1990 


28 


520 


14 


IMNor 


1920,2002 


82 


0.103 


80 


TPyx 


1890,1902,1920,1944 


12 


0.076 


62 




1966, 2011 









Notes: (a) refer ences for e ach objects are as follows. CI Aql: outburs t year dCollazzi et al. II2009I) . 
forb {M ennicken t & _Honev cutt 1995). V394 CrB: P^,^ dSchaeferllioogh . V2487 Oph: outburst year 
dPagn otta et al. 2009, who suggest the recurrence period to be possible 18 yr instead of 98 yr from dis- 
cove ry frequency), P ^^ (Schaefer 2010,). USco : Schaefer (2010) suggeste d 19 17 and 1969 eruptions, 
Porb iSchaeferl20"T(ih . T CrB: Kry, fcines. Lines & McFaul 1988). RS Oph: Isdia efer (2010) suggested 
1907 and 1 945 outburst, dSrandi et al. "200 9J V74 5 Sco: ( Schaefer "2009). V3890 Sgr: 
( ISchaeferll2 009). IM Nor: P^.b (Woudt & Wamer1l2003h . T Pyx: P„,i, ( Schaefer et al. 199Z) . 

(b) shortest recurrence period. 

(c) all ti are taken from lSchaefer ] (l2010h . 



The recurrent novae are divided into three groups depending on the type of companion star: 



1. Slighdy evolved MS stars: CI Aql, V394 CrA, V2487 Oph, and U Sco. 

2. Red-giant stars: T CrB, RS Oph, V745 Sco, and V3890 Sgr 

3. Red dwarf stars: IM Nor and T Pyx. 



Here "slightly evolved MS" means that the companion star is a main-sequence star that has 
evolved off the zero-age main-sequence (ZAMS) and its radius has expanded a little. These 
three groups a re sometimes called the U Sco, T CrB (or RS Oph), and T Pyx subclasses (e.g. 



Warner 1119951) . respectively. Table [Tjlists the recorded outburst year, shortest recurrence period. 



orbital period, and time for all the Galactic recurrent novae ever reported, where t,„ (m - 2 
or 3) is the day during which a nova decays by ;«-th magnitude from optical maximum. Table 
12] shows the WD mass estimated from the light curve analysis explained later in Section [3] f2 
time, distance d, extinction E{B - V), apparent V magnitude at maximum mv.max, and absolute V 
magnitude Mv,m-dx at maximum, and difference between the absolute V maximum observed and 
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Table 2. WD mass in recurrent novae and related objects 



object " 


WD mass 




distance 


E{B - V) " 


mv,mnx 




AMi/MMRD 




(Mo) 


(day) 


(kpc) 








(1), (2) 


CI Aql 


~ 1.2 


25.4 


4.7 


0.85 


8.8 


-7.2 


-0.3, -0.4 


V394 CrA 


> 1.37 


2.4 




0.25 


7.2 






V2487 Oph 


> 1.37 


6.2 




0.76'* 


9.5 






USco* 


> 1.37 


1.2 


6.7 


0.35'* 


7.7 


-7.5 


-3.0*, -1.5 


TCrB* 


> 1.35 


4.0 


0.96 


0.065 


2.5 


-7.6 


-2.4*, -1.3 


RS Oph 


~ 1.35 


6.8 


1.6 


0.75'* 


4.8 


-8.5 


-0.9, -0.3 


V745 Sco 


> 1.35 


6.2 




0.82'* 


9.4 






V3890 Sgr 


> 1.35 


6.4 




0.56'* 


8.1 






IM Nor 




50.0 


3.4 


0.80 


8.5 


-6.7 


-0.3, -0.4 


TPyx 




34.6 


3.2 


0.25 


6.3 


-7.0 


-0.2, -0.3 


V407 Cyg* 


> 1.37 


5.9 


2.7 


0.50 


7.1 


-6.6 


-3.0*, -2.3 


V838 Her 


~ 1.35 


1.4 


4.0 


0.43'* 


5.4 


-8.9 


-1.6,-0.1 


V2491 Cyg 


~ 1.3 


4.8 


13.0 


0.23 


7.5 


-8.8 


-0.9,-0.1 



Notes: (a) references for each objects are as follows. CI Aql: WD mass (Hachisu &Kato"'2001a'), 
distance (Hachisu & Kato 2012b). V394 CrB: WD mass ( Hachisu & Kato 2000). USco: WD mass 
dHachisu et al. 2000 aD but a su per-Chandrasekhar m ass WD is suggested by iHachisu et al. I ( l2012al) . 
distance {Ha chisu et aLll2000iJM. T CrB: WD mass iHachisu & Kato IbOO lb*). dista nce and E(B - V ) 
dBelczvnsM & Mikol aiewskall998l) . RS Oph: WD mass dHachisu et al. 120 06: Hachis u & Katol200d) . 
distance (fHiellmins et al. "1986'). V745 Sco: WD mass (Hachisu & Kato .1 2001 bl) . V3890 Sgr: WD 
mass (Hachisu & Kato 2 001b). T Pyx: dista nce ([Ha chisu & Kato 20 12^). ¥40 7 Cyg: W D mass 
dHachisu & KatoH2012al). and dista nce (Mun ari et al. |[2012h. E(B - V) (Shore et aLll201ll). V838 



culated from the results of 



Her: WD mass ( Kato. Hac hisu & Cassatellai200^ ). h ( Harrison & Stringfellow ul994j) . distance is cal- 



iKato et al. I d2009l) and th e present value of E(B - V). V2491 Cyg: WD 
mass (Hachisu & Kato 2009), (Munari et al. 2011'), distance is calculated from the same method in 
lllachisu & Kato ( 2012b), E(B - V) (Munari^ et al. 20 11). 

(b) all t2 for recurrent novae are taken from lSchaeferlfcOlOl) unless otherwise specified. 

(c) all E(B - V) for recurrent novae are taken from lSchaefer1(l2010l ) unless otherwise specified. 

(d) E(B - V) is obtained from the Galactic dust absorption map of NASA/IPAC (http:// 
irsa.ipac.caltech.edu/ applications/DUST/). 



the V maximum calculated from the MMRD relation. We added three related objects in order to 
compare them with the Galactic recurrent novae. Recurrent novae in the LMC and M31 are not 
listed in these Tables. 

Note that in order to classify a nova into the subgroup of recurrent novae, just a rapid de- 
cline in the light curve, which suggests a massive WD, is not enough. For example, the clas- 
sical nova V838 Her showed a light curve as steep as that of U Sco, also the light curve of 
V2491 Cyg decayed similarly t o that of RS Oph. Their WD masses are esti mated to be very 
massive, 1 .35 Mq for V838 Her jKato et al. 12009) . and 1 .3 Mq for V2491 Cyg dHachisu & Kato ' 
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12009"). However, these nova ejecta showed heavy elem ent enrichment, Z = 0.12 for V838 her 
IVanland ingham, Starrfield & Shore 1997; Sc hwarz et al. 2007) and Z = 0.14 for V2491 Cyg 
dMunari et al. I l201ll) . suggesting a long recurrence period between nova outbursts during which 
diffusion processes work to dredge up WD core material into hydrogen-rich envelope. V2491 
Cyg also showed a much shorter SSS phase (10 days) than that of RS O ph (60 day s). A relat- 
ively long duration of a SSS phase stems from a mass-increasing WD (iKatoll20li see Section 
13.31 ) while a relatively short one indicates a mass-decreasing WD as in normal classical novae. 
Therefore, both V838 Her and V2491 Cyg are highly likely classical novae rather than recurrent 
novae. 



3. Theory of Nova outbursts 
3.1 Stability of shell flashes 

A mass-accreting WD has a geometrically thin envelope before a nova eruption. The envelope is 
usually cold because of radiative cooling. As the envelope mass increases with time, the bottom 
of the envelope is gradually heated up due to compressional heating (gravitational energy release). 
When the envelope mass reaches a critical value, the radiative cooling almost balances with the 
total release of gravitational and nuclear burning energies. In such a situation, even a very small 
increase in the temperature causes a large increase in the nuclear burning energy release. This 
increase in the energy release cannot be dissipated away by radiation, and so results in a runaway 
of hydrogen nuclear burning. This is the beginning of a nova outburst. 

The ignition mass of hydrogen-rich envelope is obtained from the condition that the radiative 
cooling is balances with the total release of gravitational and nuclear burning energies. This 
ignition mass depends on the WD mass, mass accretion rate, WD temperature, and opacity (i.e., 
chemical composition) of the envelope. In recurrent novae, the ignition mass does not depend 
on the thermal history of the WD because the recurrence periods are shorter than ~ 100 yr The 
ignition mass is smaller for a more massive WD because of its large compressional heating. Fig. |2] 
schematically depicts ignition masses by solid lines in the region of Mace < M^t, where hydrogen 
shell burning is unstable. 

During nova outbursts, all of the accreted matter will be blown off. Moreover, a part of the 
WD core is eroded due to diffusion of hydrogen into the core (see the next subsection). Then, the 
WD mass decreases after every nova outburst. In the case of recurrent novae, however, hydrogen 
shell flash is relatively weak and a certain amount of processed helium is added to the WD after 
every outburst. Therefore the WD mass increases (see Section l53b . This recurrent nova region is 
indicated by shadow in Fig. |2] Their accretion rate is as large as Mace ~ 10"^-10~^Mo yr~' and 
the WD is quite massive (Mwd S I.IMq). 

When the mass accretion rate is larger than Mace ^ Mjt, hydrogen nuclear burning is stable 
and ash of hydrogen burning accumulates on the WD. When the mass accretion rate exceeds Mcr, 
the envelope expands to blow optically thick winds. The system undergoes an accretion wind 
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Figure 2. An evolutionary path of typical SN la progenitors (red solid line with arrows) on the map of 
WD responses to the mass accretion rate M.^^c- The progenitor evolves from the optically thick wind phase, 
through SSS phase (steady H-buming phase) and recurrent nova phase (H-shell flash phase), and finally 
explodes as an SN la at star mark. In the region above the line of (Mace > ^cr)> strong optically thick 
winds blow, which stabilize binary evolution in the SD scenario. On the other hand, a common envelope 
is formed in the original DD scenario, which results in a binaiy consisting of a double degenerate (WD) 
system. In the region of Mjt ^ Mice ^ Ma, we have steady hydrogen shell burning with no optically thick 
winds. There is no steady-state burning below the line of M,, (M.^^^ < Mj,). Instead, hydrogen shell-burning 
intermittently occurs and results in a nova outburst. The envelope mass AMi„, at which an intermittent 
hydrogen shell flash ignites, is also shown. Shadow area schematically shows the region in which the 
recurrence period is shorter than 100 yr, that is, the region of recurrent novae. The original data are taken 
from lNomotol lll982h . lHachisu & Katol j200ld) . and lHachisu et al.l feOIOh . 



evolution as illustrated in Fig. [T]while the binary was supposed to undergo a common envelope 
evolution in the original DD scenario. 



Stability analy sis of hydrogen shell - burnin g had been established a long time ago based on 



1983). His results were confirmed recently with the 



These sta bility criteria were also confirmed directly i n 



linear analyses bv ISienkiewicz I (11975 
OPAL opacity (Nomotoetal T l2007h 
many numerica l calcii lations of shell flashes (Prialnik & Kovetz II1995I : iJose & Hernanz 1119981) . 
Starrfield et al. 1(120041) . however, presented an opposite result that all accreting WDs are thermally 
stable (this means that no novae/recurrent n ovae occur), thus , that the accreting WDs always grow 



in mass to reach the Chandras ekhar mass (IStarrfield et al. 112004) . Their calculations were criti 



cized bv iNomoto et al.l (l2007h for the reason that their results are artifacts originatin g from the 
lack of sufficient numerical grids. Recently, the same group ("Starr field et al. presented 
a completely different result to their previous one that all accreting WDs are unstable for shell 
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flashes and virtually no SSS phase exists in classical novae. Again their results are not sup- 
ported by the stability analysis cited above as well as many observed SSS phases of classical 
novae/recurrent novae (see, e.g. jHenze et al. boioiboiil for an M3 1 nova survey). 



For the last, we must note an inappropriate statemen t on the "proper pre ssure," which is 



widely used in estimating the ignition mass of novae (e.g. lOsborne et al. 11201 11) : there is a sup 



posed simple relation between the pressure at the bottom of an envelope and the ignition mass, 

cm , where P* is called the "proper pressure," 
MwD is the WD mass, Me,„ the hydrogen-rich envelope mass, /?wd the WD radius. Thus the 
ignition mass is in dependent of the mass-accretion rate. This relation may originate from a form- 
alism proposed by Shara ( 198ll) and the numerically obtained critical pres sures which were b ased 
on a simple env elope model of polytropic/adiabatic approximations bv iMacDonald I (119831) and 
Fuiimotd (119821) . This simple expression, however, has discrepancy of envelope mass as large as 
an order of magnitude or more. Numerical results obtained with non-adiabatic, dynamical calcu- 
lations with the OPAL opacity clearly show that the ignition con dition is largely different from 
the "proper pressure" and depends on the mass accretion rate (e.g. Prialnik & Kovetz1ll995 ). 



3.2 Chemical composition of ejecta 

Classical novae show heavy element enrichmen t in their ejecta by an amount of a few to several 



tens percents by mass (e.g., iGehrz et al. II1998I) . The enrichment is interpreted as dredge-up of 



WD core material into ejecta: Hydrogen diffuses into WD cores during the accreting phase and 
ignites below the original WD core surfaces, thus a part of core material is mixed by convection 
into the hydrogen-rich envelope and blown off in ejecta d Prialnik Ill986l) . The diffusion occurs 
when Mace 5 10"^ Mq yr"', i.e., in a long timescale of Menv/M^cc ~ lO""* Mq/M^^c <: 10^ yr for 
classical novae. Therefore, the diffusion process is not effective in recurrent novae, in which the 
recurrence periods are as short as 10-100 yr Thus, heavy element enhancement is not expected. 

Recurrent novae, however, possibly show hydrogen depletion. In the very beginning phase of 
a nova outburst, convection widely develops and carries processed helium up into the envelope. 
This convective mixing reduces the hydrogen content by 10%-20% by mass for massive WDs 
like in recurrent novae. Therefore, we expect the hydrogen content of X x 0.50 and the solar 
metallicity of Z « 0.02 in envelopes of recurrent novae. 

Because dredge-up of WD core material is not theoretically expected in recurrent novae, it 
is hard to know whether WDs in recurrent novae are made of carbon and oxygen (CO) or of 
oxygen, neon, and magnesium (ONeMg). This is important from the SN la progenitor point 
of view, becau se only CO WDs can exp lode as an SN la while ONeMg WDs will collapse to 
a neutron star ( Nomoto & Kondj 1991 ). In the SD scenario, recur rent novae con tain either a 



CO WD (increased from a initial WD mass of < 1.07 Mq, see e.g.. lUmeda et al. 1999^ for an 



upper limit of CO WD masses) or an ONeMg WD. In the DD scenario, all of the massive WDs 
(> 1.07 Mq) would be an ONeMg WD, because WDs hardly grow in mass during and after a 
common envelope phase. 
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Recently, Mason ( 201 claimed that U Sco harbors an ONeMg WD because of a high Ne/O 
line ratio in spectra of the 2010 outburst, comparable to those in several neon novae, but much 
larger than those in CO novae. Unfortunately there are no such data in other recurrent novae 
other than classical novae. As we show later in Section |4] recurrent novae have a very different 
evolutional history from those of classical novae. Mass-increasing WDs like in recurrent novae 
develop a thin helium layer underneath the hydrogen burning z one. After the helium sh ell grows 
in mass to reach a critical value, a helium shell flash occurs (e.g..'Kato & Hachisu 1 1981 . Helium 
burning produces a substantial amount of Ne and Mg (Shara & Prialnik I I 19941) . Even if the WD 
is made of carbon and oxygen, the hydrogen-rich envelope could be contaminated with neon 
and magnesium. Therefore, some recurrent novae would show neon enhancement if it occurs 
shortly after the latest helium shell flash, because helium accumulation after every recurrent nova 
outburst could hide or dilute the products of helium shell flashes. Large Ne/O line ratios of ejecta 
do not always indicate an ONeMg WD. Instead, it may indicate that a CO WD develops a helium 
layer during e very recurrent nova outburst. We suggest that the high Ne/O line ratio in U Sco 
( Mason |[2011 ) is also an evidence of mass-increasing WD, although we need more observational 
information on the chemical compositions of U Sco and the other recurrent novae. 



3.3 Light curve analysis of recurrent novae: how to determine white dwarf mass 

After hydrogen burning sets in, the envelope rapidly changes its structure from a geometrically 
thin shell to a bloated spherical configuration. As the envelope expands from the WD surface 
to a red-giant size, the photospheric temperature rapidly decreases down to log T (K)~ 4.0 or 
lower. This spherically extended configuration is stable against hydrogen shell burning. Then the 
nuclear burning settles in a steady-state. The envelope mass decreases owing to nuclear burn- 
ing and wind mass-loss. The hydrogen burning continues until the envelope mass decreases to 
reach a critical value for the extinguish condition. In this decay phase, the photospheric radius 
gradually shrinks whereas matter goe s away. So the photospheric temperature increases with time 



dKato&Hachisull 19941 : iKatoll 19991) . After the wind stops, the envelope mass still continues to 



decrease owing to nuclear burning. This stage corresponds to a SSS phase because the temperat- 
ure is high enough to emit supersoft X-rays. At the termination of nuclear burning, the envelope 
structure comes back to a geometrically thin shell. 



In the decay pha se of nova outbursts, fre e -free emission of optica lly thin ejecta dominates 
continuum flux (e.g.. lGallagher & Nev iHachisu & Kato I ( 12006 ) model ed free-free emis 



sion li ght curves of novae based on the optically thick wind model calculated bv^ Kato & Hachisu 
(I1994I) . The decline rate of the model light cur ves depends sensitively on the WD mass and 



weakly on the chemical composition of envelope (IHachisu & Kato Il2006h . Especially, when the 
WD mass is close to the Chandrasekharmass limit Mch ~ 1 .4 Mq, the Ught curve depends sharply 
on the WD mass. This is because the WD radius is very sensitive to the increase in mass near 
Mch- Therefore, we can in principle determine the WD mass by fitting the model light curve with 
the observed one (see, e.g. . iHachisu & Kato 



Here, we should note that dynamical calculations of nova outbursts have numerical diffi- 
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culties, especially in the treatment of surface boundary condition. Nariai et al (1980) already 
pointed out that expansion velocities are very different among different groups (see their Table 
4). We note that Israel code ( Idan. Shaviv & Shavivll2012h tends to derive violent outbursts with 
large expansion velocities and large mass-loss rates, due to inadequate boundary conditions in- 
stead of solving the envelope up to the photosphere. No dynamical code has been succeeded so 
far in calculating extended stages of nova outbursts and their light curves. On the other hand, 
the optically thick wind theory was developed in order to calculate such extended stage of nova 



outburst and has been succeeded in reproducmg a number or nova light cu rves (Hachis u et al. 
20061; iHachisu & Kato|[200^l2007l : [Hachisu. Kato & C assatella' 2008[|Kato et al. II2009I) as well 
as expansion velocities (see Figure 18 in lHachisu & Katoii2007l) . 



One of the characteristic properties of recurrent nova light curves is a mid-plateau phase, 
that lasts 10-60 days after the brightness decayed by several magnitude. The plateau phase often 
overlaps with a SSS phase, i . e., it begins/ends when the supersoft X -ray flux rises/decays (e.g., 
Hachisu. Kato & Lunal 20071 : Osborne et al. 2011 : Ness et al. 20121) . S uch plateau phases can 



be explained by the contribution of an irradiated disk as demonstrated bv IHachisu et al. I (l2000ah 
for U Sco and by Hachisu et al. (2006, 2007) for RS Oph. The presence of an accretion disk is 
also suggested in T CrB, in which the secondary max imum in its light curve can be interpreted as 
the contribution from an irradiated large tilting-disk (IHachisu & Kato On the other hand, 

in classical novae, such mid-plateau phases are rarely reported. The contribution of disks is also 
suggested in classical novae, but it appears only in a very late phase of nova outbursts (see, e.g., 
Hachisu. Ka to & Katd'2004l, for V1494 Aql). It is because the orbital period is usually much 
shorter in classical novae than in recurrent novae, and thus classical novae have much smaller 
disks that do not contribute much to optical light curves. 



We think that accretion disks are not entirely blown off during the outburst in recurrent novae. 
One of the indications is the presence of these mid-plateau phases. All the recurrent novae shows 
more or less a plateau ph ase. The other indications are the detection of flickering 8 days after 
the U Sco 2010 eruption (IWorters et al. i2010) and formation of accretion disk at least 8 days 
after the 20 1 outburst of U Sco ( Mason et al. 112012 ) . During the mid-plat eau phase, orbital ligh t 
curves of U Sco have asymmetric multiple-peak shapes outside the eclipse (Schaefer et al. 20121), 
which can be reproduced by spiral arm structure on the accretion disk (see, e.g., H achisu et al. I 
2004 for V1494 Aql). Therefore, it is very likely that the disk was not gone but survived during 
the 2010 outburst. On the other hand, some numerical cal culations i ndicated that the accretion 
disk is entirely disrupted by the fast moving ejecta (D rake & Orlando 2010). This seems to 
be inconsistent with the presence of the disk only 8 days after the outburst. We suppose that 
numerical calculation cannot resolve the accretion disk that has high density region near the 
equator. 



Outburst light curves of recurrent novae are a summation of each contribution from a WD, 
irradiated a ccretion disk, and irrad iated com panion. Such composite light curve s are modeled for 
RS Op h bv lHachisu et al. I (l2006h. U Sco bv IHachisu et al. I (l20()0ah. T CrB bv IHachisu & Kato 
(ll999), V394 CrA bv IHachisu & Katol ll200nh . and CI Aal bv lHachisu & Kato I l l2001a '). The 
resultant WD masses are listed in Table |2] 
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3.4 Maximum magnitude versus rate of decline relation 

There is a statistical relation between the maximum V magnitude, My rnax, and the rate of decline , 
f2 or f3, in optical light curves of classical novae (e.g.. lSchmidtll957HDella Valle & Liviolll995h . 
This is called the MMRD relation. Distribution of individual novae in the /3-Mv;max plane shows 
a main trend fro m the left-top to the right-bottom but there is a lar ge scatter above/be l ow the 



main trend (e.g., lOownes & Duerbeck I Hoool: iKasUwal et 31711201 ih . iHachisu & Kato I ( l2010h 



theoretically explained the main trend and scatter of the MMRD relation. The main parameter 
that determines the main trend is the WD mass, and the second parameter for the scatter is the 
ignition (initial envelope) mass, in other wards, the mass accretion rate to the WD. If the mass- 
accretion rate to the WD is relatively larger, the ignition mass is smaller (see Fig. |2]l, so that the 
peak brightness is fainter. This second parameter c an reasonably explain the scatter of individual 
novae (see, e.g.. Fig. 15 of lHachisu & Katoll201(]l) . 

It is, however, frequently discussed that the MMRD relation cannot be applied to recurrent 
novae. Table |2] show the difference between the absolute V magnitude and what the MMRD 
predicts, that is, AMvmmrd = Mvmmrd - A^v.max, for six recurrent novae and three very fast 
classical novae. Here we estimate the absolute V magnit ude Mt^max using the d i stance and ex- 
tinction in Table m For the MMRD relation, we use (1) 'Dow nes & Duerbeckl (|2000|) . that is, 
Mvmmrd - -10.79 + 0.92 -i- (1.53 ± 1.15)logf2 for faster novae (logf2 < 1.2), but Mv^mmrd - 
- 8.71 + 0.82 + (1.03 + 0.5 1) log fa for slower novae (log t2 > 1 .2), and also, for comparison, (2) 



Delia Valle & Livio 



1995 '). that is, Mv.mmrd = -7.92 - 0.81 arctan((1.32 - log f2)/0.23). Four 
recurrent novae, CI Aql, RS Oph, IM Nor, and T Pyx show reasonable agreement (within 1 cr 
error) with what two MMRD relations (1) and (2) predict. Two classical novae V838 Her and 
V2491 Cyg also show good agreement with what MMRD relation (2) predicts. 



On the other hand, U Sco and T CrB have a very similar absolute V magnitude at maximum, 
Mv,max = -7.5 and -7.6, respectively, and both of them show a large deviation from the MMRD 
prediction. Even if we use MMRD relation (2), the absolute magnitude of U Sco is still 1 .5 mag 
fainter than the prediction. The symbiotic classical nova V407 Cyg also shows a very large devi- 
ation from the MMRD relation. We think that these large deviations from the MMRD relation are 
one of the indications of super-Chandrasekhar mass WDs (Section lCTT ) because a large deviation 
means a small envelope mass at the optical peak (i.e., small ignition mass), which is a strong 
indication of very massive WDs as explained in section IXTI 



4. Evolutionary status of recurrent novae 



In this section, we discuss the evolutionary status of recurrent nova systems and how they will 
evolve toward SN la explosions. In the original DD scenario, there are no known paths to recur- 
rent novae especially for the RS Oph type systems. In the SD scenario, there are two channels to 
SNe la, i.e., the WDh-MS (or simply MS) channel and WD+RG (or symbiotic) channel. In the 
MS channel, a pair of MS stars evolve to a binary consisting of a helium star plus an MS star after 
the first common envelope evolution (e.g., Hachisu et al.ll999b ). The primary helium star further 
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Figure 3. Initial and final binary parameter ranges proposed by a recent SD scenario jHachisu et alj 
l2012bl) . Here, Mo is the companion mass and the orbital period. An initial binary system inside the 
region encircled by a solid line (labeled "initial") is increasing its WD mass up to the mass of Mwd = 1-38, 
1.5, 1.6, 2.1, and 2.2 Mq (from outside to inside) and then reaches the regions labeled "final" when the 
WD stops growing in mass. Currently known positions of recurrent novae and supersoft X-ray sources are 
indicated by star mark (★) for U Sco (e.g.. .Hachisu et aO2000ah . filled square for V Sge jHachisu & Kato I 
l2003bh, filled triangle fo r T CrB (e.g.. lBelczvnski & MikolaiewskalFl993) . and open diamond for RS Oph 
(e.g.. lBrandi et aL l2009l) . In the very last stage, the companion mass is reduced to 1—2 Mq for the WD+MS 
systems, and < 1 for the WD+RG systems. 

evolves to a helium RG with a CO core and fills its Roche lobe followed by a stable mass-transfer 
from the primary helium RG to the secondary MS star The secondary MS star increases its mass 
and is contaminated by the primary's nuclear burning products. The primary becomes a CO WD 
after all the helium envelope is transferred to the secondary MS star. Then the secondary evolves 
to fill its Roche lobe and mass transfer begins from the secondary MS (or subgiant) star to the 
primary CO WD. When the mass-transfer rate exceeds M^-, an accretion wind evolution is real- 
ized. The mass-transfer rate gradually decreases to below M^i and the optically thick winds stop. 
Then all the accreted matter is burned on the primary CO WD. The binary enters a persistent 
supersoft X-ray source (SSS) phase. We should note that we must distinguish this persistent SSS 
phase from transient SSS phases in classical novae and recurrent novae. The mass-transfer rate 
further decreases to below M^t and the binary enters a recurrent nova phase. 

In Fig. [3] binaries within the area labelled "initial" evolve downward and reaches the "fi- 
nal" area when its mass increases to Mwd > 1 -38 Mq. The supersoft X-ray source V Sge and RX 
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J05 13 . 9-695 1 are regarded as a binary corresponding to the accretion wind phase (iHachisu & Kato 
2003allbl) . i.e., these binaries are now evolving toward SNe la. The present status of U Sco and 
other recurrent novae in Table [1] corresponds to the final stage of the above three phases (wind, 
SSS, and recurrent nova). If a WD explodes during the wind phase, the resultant SN la shows 
some indications of the interaction between ejecta and circumstellar matter (CSM) and the pres- 
ence of hydrogen, because the WD winds and matter stripped from the companion should form a 
rather massive CSM. If it explode during the SSS or recurrent nova phase, there are weak indic- 
ations of hydrogen because the secondary MS mass are ~ 1-2 M© and their size is compact. In 
multidimensional dynami cal calculatio ns of SN la explosion, however, a very massive compan- 
ion is often assumed (e.g. lKasen" 20ld 2-6 Mq MS) and they conclude that the SD scenario is 



inconsistent with no UV excesses as well as no detections of hydrogen lines. However, we think 
that these secondary masses are rather large and not realistic. 



In recent SN la pr ogenitor models of lJustham I (1201 11) . iDi Stefano. Voss & Claevi (1201 lb . and 



Hachisu et al. (1201 2bl) . they suppose that a rapidly rotating WD does not explodes as an SN la 



even if it exceeds Mwd = 1-38 Mq, because the central density is too low to ignite carbon. They 
further assume that it takes a long time before the SN la explosion. Once it exceeds Mwd - 
1.38 Mq, then the explosion is postponed until it spins down and the central density increases 
to reach a critical density for carbon ignition. During such a long waiting time, the companion 
evolves off the main-sequence and become a WD or transfers its mass to the primary WD and 
becomes a cataclysmic variable with a short orbital period. When the primary WD explodes as 
an SN la, neither hydrogen nor a bright e x-companion is left, but only a rather dark WD or a red 



dwarf ex-companion is left. IHachisu et al . (2012b) showed that such a dark ex-companion is in a 



majority of SN la remnants coming from the MS channel. 

In the symbiotic channel we start from a very wide binary in which the primary evolves to 
an asymptotic giant branch (AGB) star. The binaries unde rgo a common envelop e-like evolution 



during the superwind phase of the primary AGB star (e.g.. IHachisu et al. Ill999al) . The separation 
of the binary shrinks to the orbital period of 30 - 800 days and it becomes a pair of a CO WD and 
an MS star. Then the secondary MS evolves to an RG with a helium core and fills its Roche lobe. 
Then the system enters an accretion wind evolution. After that, the binary becomes a persistent 
SSS, and evolves finally to a recurrent nova. Depending on the binary parameters, the WD mass 
reaches Mwd = 1-38 Mq during one of the three phases (wind, SSS, and recurrent nova). One 
of such evolutionary paths is plotted in Fig. |2]by a red solid line. The SMC symbiotic X-ray 
source SMC3 corresponds to the SSS phase and symbiotic recurrent novae like T CrB do the 
last phase. If the WD explodes as an SN la during the first two phases, hydrogen should be 
detected in th e SN la explosion, because substantial amount of matter is stripped out from the 



RG (see, e.g.. lKato. Ha chisu & Mikolajewska 2012b,, for SMC3). In the recurrent nova phase. 



however, the companion has already substantially exhausted its envelope through the long-lasted 
mass-transfer phase. Its envelope mass may be small (< 0.5 Mq), because they have developed a 
hehum core. 



In the spin-up/spin-down scenario of SNe la, the SN explosion could be postponed until the 
secondary RG evolved off the red-giant branch and became a helium (or CO) WD. These dark 
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companion WDs could not be detected just before and after the SN la explosion. iHachisu et al 
(|l)12b) showed that in their rapidly spinning WD progenitor models, these dark companions are 
the majority of SN la progenitors. Therefore, we cannot expect any indication of hydrogen or 
bright companions in such SN la explosions. 



Observations of SNe la have provided the following constraints on the nature of compan- 
ion st ars. Some evidences support the SD model, such as the presence of circumstellar mat- 



ter jPatat et al. Il2007t [Sternberg et al. 11201 ll Foley et al. 1120121) and detections o f hydrogen in 



the ci rcumstellar matter-ejecta in teraction type SNe (type la/IIn) l ike SN 2002ic (Ham uy et al 



20031) . SN 1604 (Kepler's SNR) dChiotellis. Schure & Vinkl(20 1 2 l). and PTFllkx ( Dildav et al 
20121) . On the other hand, there has been reported no direct indication of the presence of com- 



panions, which is usually thought to be unfav orable for the S D scenario, e.g., (1) the lack of 
companion stars in th e ima ges of SN 20 life ( Li et al. 201 lb. some SN la rem nants (SNRs) 
(ISchaefer & PagnottaJJ^l2), SN 1572 (Tycho's SNR) jKerzendorf e t al. ||2009|) and SN 10() 6 
(iKerzendorf et al. Il2012h . (2) the lack of UV excesse s of early-time l ight curves ( iKasenlEoioh . 
and (3) the lack of hydrogen features in the spectra (iLeonari 2007). Both (2) an d (3) are ex- 
pected from the collision between ejecta and a companion. iHachisu et al. I (l2012bl) showed that 
their SD scenario mainly produces unseen companions at SN la explosion because of a long 
spin-down time. They also showed that the WD mass distribution at SN la explosion predicted 
by their model agrees well with the distribution of SN la brightness observed if the brightness is 
correlated to the WD mass. These results strongly support that a major route to SNe la is the SD 
system and they explode through a recurrent nova stage. 



5. Elementary processes in binary evolution 



The SD scenario is based on different elementary processes from those assumed in the original 
DD scenario, which result in a very different binary evolution as well as population synthesis. 
Thus, it is important to clarify theoretical backgrounds of each elementary process. In this section, 
we critically examine each elementary process based on the results of nova studies, because many 
of these ingredients are closely related to nova phenomena. 



5.1 Common Envelope Evolution 

Common envelope evolution is a fundamental process of binary evolution for both the SD and 
DD scenarios. The first common envelope evolution in a binary occurs when the primary star 
evolves to a red-giant and fills its Roche lobe. Common envelope evolutions were formulated 
in a very simple form based on the total energy conservation (gravitational energy of the orbit 
plus envelope binding energy) before and after the common envelope evolution. The ratio of the 
binary orbit before (init ial) and after (final), Af/Aj, is given by a simple form using an efficiency 
parameter a(< 1) (e.g. Ilben & Livio I ll993h . A smaller value of a leads to a smaller ratio of 



Af/Ai, i.e., a larger shrinkage of the binary orbit. In the original a-formalism, a is only the para- 
meter Some authors had a different expression to calculate the gravitational energy of the RG 
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envelope by using /l-parameter (IWebbinklll984t Ide Kool I ll990h . which represents the effective 
mass-weighted radius of RG stars . Then the product of aA behaves like the original a. Obser- 
vationally, IZorotovic et al. I (1201 Ol) derived a — 0.2 -0.3, assuming that a is universal together 



with A - 0.5, while Rebassa-Mansergas et al. ( 201 2h obtained a ~ 0.25. iDe Marco et al. I (1201 ll) 
obtained the dependence of a on the mass ratio q - M2/M1 such as a ~ 1 for q - 0.1, a ~ 0.2 
for q - 0.3, and a < 0.1 for q = I. These new a-parameters work together with the effective 
radius par ameter of A, for exampl e, A - 0.4 for Mi - 5 Mq or A - 0.45 for Mi -1 Mq of AGB 



stars (e.g.. iDe Marco et al. 11201 11) . These recent analyses suggest that aA ~ 0.1-0.2. In many 
population synthesis calculations that support the DD scenario, however, a larger value of aA ~ I 
was preferentially used independent ly of q, b ecause such a larger value yields a larger birth rate 
of SNe la in their DD models (e.g.. lRuiter e t al. 2009). It is however interesting that the birth 
rate of SNe la coming from SD binaries increases and becor nes comparable to t he rate coming 
from DD binaries if we have a smaller value of aA < 0.5 (e.g..^ Ruiter et al. 112009 ^. This suggests 
that if we adopt a realistic value of a A such as a A ~ 0.1, the SD binaries dominate the progenitors 
of SNe la in population synthesis. 

One of the recent t rends in the population sy r ithesis community is to assume a n extremely 
large value of a (> 1) dHurlev. Tout & Polsll200a iTout 1 120051; iMen^ & Yang1l2012h . If there is 
a large additional energy source (e.g., recombination energy), the binary separation will even in- 
crease after the common envelope evolution (IWebb ink 2008 ). However, such a mechanism seems 
to be very difficult to work effectively because it needs all the envelope matter instantaneously 
recombined in the accelerating zone and the released energy efficiently conve rted to the binding 
energy without radiative loss (see, e.g., a criticism bv lSoker & Harpaz 1 l2003h . This id ea may be 
introd uced in order to reproduce symbiotic recurrent nova systems in the DD scenario dWebbink 
20081). However, recent observational estimates, a A ~ (j.l for AGB stars, do not support such an 
extremely high value of aA (> 1) ("e.g.. IZorotovic et al. l2010h . 



Nelemans et al. I (l2000l) proposed a new formalism for common envelope evolution based on 



the assumption that angular momentum loss is proportional to the amount of lost mass, AJ/ J 
-yAM/M. This formalism c ontains a p roportionality parameter y, so dubbed "y- formalism." This 



rop' 



forma lism is criticized by IWebbink I (12008) as it does not actually limit Af. IZorotovic et al. 



(l201flh also criticized the y-formalism because the a-formalism is much more reasonable from 
the statistics of a large number of binary parameters. 



5.2 Accretion wind evolution versus common envelope evolution 

After the first common envelope evolution, the binary becomes a pair of a WD and a secondary 
star After some time elapses, the secondary evolves to fill its Roche lobe. Then the mass transfer 
begins from the secondary to the primary WD through the Ll point. Infalling matter forms 
an accretion disk and finally accretes on to the WD. The accreted matter quickly spread over 
the entire WD surface and hydrogen nuclear burning ignites after some critical amount of mass 
accumulates. If the secondary is an evolved RG more massive than the WD, the mass accretion 
rate exceeds the critical value, M^r (see Fig. |2]), and strong optically thick winds blow. Once the 
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wind occurs, the wind carries away mass and angular momentum. Then the orbital separation and 
the size of the secondary's Roche lobe change in time. The mass-transfer rate is regulated to keep 
the secondary's radius equal to the Roche lobe radius. If t he mass-transfer rate is smaller than 



Mace ^10 Mq yr , a common envelope does not form (I Hachisu et al. Ill999at iHachisu et al 



[i999b). This is the accretion wind evolution (Fig. [Til. When the accretion rate is very large 
A^acc ^ 10"^ Mq yr however, the photospheric radius of the WD expands over the Roche lobe, 
and a common envelope evolution may occur. 

If the accretion wind evolution is taken into account, the condition for the second common 
envelope evolution is relaxed to q < 1.15 from q < 0.79, where q - Mrg/Mwd- The former 
condition is obtained with the mass and angular momentum conservations including those of the 
winds, and the latter without the winds. Therefore, some binaries can avoid a formation of the 
second common envelope and keep the binary separation. The binary still consists of a mass- 
accreting WD and a lobe-filhng secondary. The WD can grow in mass during the accretion wind 
evolution dHachisu et al. Ill996h . 

This is the first step that widens the binary parameter range of SN la progenitors. It is 
interesting to note that optically thick winds have been established in nova outbursts. The U 
Sco 2010 outburst suggested that an accretion disk formed (or still existed) in a ve ry early phase. 



only 8 days after the optical maximum (see Section [331 and iMason et al. II2012I) . In this early 



phase, the optically thick wind had not yet stopped (it ends at aro und Day 13) but th e companion 



and accretion disk emerged from the extended envelope at Day 8 (IHachisu & Kato 112012a ). This 
means that accretion from the companion can be observed when the companion emerged from 
the WD photosphere, while the optically thick wind still continues. We think that this is an 
observational example of accretion wind evolution. 



5.3 Mass stripping in accretion wind evolution 



When the accretion wind occurs as illustrated in Fig. [T] the fast wind (> 1000 km s ') hits the 
companion's surface that fills the Roche lobe. A very surface layer of the companion's hemi- 
sphere may be strippe d by the winds, which reduces the mass outflow rate from the secondary. 



Hachisu et al. I (Il999ah first incorporated this effect into their binary evolution and established the 
symbiotic channel to SNe la. This is the second step that widens the parameter range of binary 
evolution toward SNe la. 



The accretion wind evolution including the mass-stripping effect is also essential to un- 
derstand the m echanism of quas i -period ic behavior of the two SSS binaries, V Sge and RX 
J0513.9-6951. IHachisu & Katol ( l2003allbh explained both the optical and supersoft X-ray light 
curves of these two objects based on the accretion wind evolution that leads to a limit cycle of the 
light curve behavior They obtained the WD masses of ~ 1.2-1.3 Mq, and the mass-increasing 
rates of the WDs ~ 10 *' Mq yr ' for the both objects from light curve fitting. Therefore, these two 
systems are also candidates of SN l a progenitors. Ther e are several similar objects in our Ga laxy, 
the so-called V Sge-type stars (.Steiner & Diaz Ill998h . iKafka. Anderson. HonevcuttI (l2008h and 
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Linnell et alJoOOSh found that QU Car is a very similar object to V Sge and belongs to the V Sge- 



type stars. This also suggests that many such objects could be misclassified to other categories of 
variable stars and the number of V Sge-type stars would significantly increase. 



The successful results of limit-cycle light curves strongly suggest that (1) both the processes 
of "accretion wind" and "mass stripping" can be realized in binary evolution and (2) we expect 
a wider parameter region of binary evolution toward SNe la by avoiding the second common 
envelope evolution. 



5.4 AGB super winds in wide binaries 



Stars with a mass less massive than M < 8 Mq lose most of its envelope mass in a short 
timesca le (< 10^ yr) at the final stage (at its asymptotic giant branch, AGB, phase) of evol- 



ution. 



Hachisu et al. (fl999a) incorporated this AGB superwind as an important process in 



binary evolution. Due to large mass-loss rates o f superwinds (up to several xlO Mq yr 
[Groeneweg en et aL 2002 : Guandalini et al. even very wide binaries (a, = 1500 - 30000 

Rq) experience a decay of the orbital separation that is similar to a common envelope evolution. 
After that, it becomes a wide binary corresponding to the orb ital periods of symbiot ic stars. This 
is the symbiotic channel of the SD scenario first proposed bv lHachisu et al. I (Il999ah . 



The original DD scenario did not include the evolution driven by superwinds in very wide 
binaries. Many population synthesis calculations still started from zero-age binaries excluding 
very wide binaries. Thus, after the first common envelope evolution, the separation shrinks from 
fl, ~ 700-2200 Rq to af ~ 5-17 Rq {Poih ~ 0.5-1.5 days) if we use aA ~ Q.l from the recent 
trend of small values of a (see Section l5Tt . Then the resultant binary periods lie in the region 
of 0.5-1.5 days, which is not consistent with the present states of symbiotic recurrent novae but 
consistent with the MS channel of SD scenario. This is the reason why the SN la rate coming 
from SD binaries becomes large and exceeds the rate coming from DD binaries if we assume a 
small value of a A.. 



5.5 Mass increasing rate of white dwarfs 



The original DD scenario was based on the mass increasing efficiency of WDs being very different 
from that of the SD scenario. In the original DD scenario, mass-accreting WDs can grow in mass 
only if it stays in the narrow strip region (Mjt < M^cc ^ Ma) in Fig. |2] These WDs, however, 
suffer helium shell flashes, which take away a large part of the helium envelope mass. Thus the 
total efficien cy of mass-increase is very s mall. Such a small mass-inc reasing efficiency is based on 



the claim bv lCassisi. Iben & Tornambd (11998.) and iPiersanti et al. I (11999. .2000) who calculated 



shell flashes on low-mass WDs (< O.SMq) using a spherical symmetric hydrostatic code with 
the Los Alamos opacity. They concluded that the frictional process due to the companion's 
motion should be very effective to take away most of the envelope mass, simply because the 
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envelope photosphere expands beyond the binary orbit, although their calculation did not include 
any frictional effects in their codes. 



It is worth noting that ICassisi et al. (Il998 l) emphasized 'the Los Alamos opacities are very 



similar to the OPAL opacities' (see the last sentence of Section 4 in their paper). However, this is 
clearly an incorrect statement ( see a comparison betwe en nova envelopes calculated with the new 



and old opacities in Fig. 15 of Kato & Hachisu II 19941) . As we have already known, calculations 
with the old opacity could derive very different conclusions. Even in the sa me framework of 



1980's phys ics, h o wever, th eir claim was not supported by the results given bv lKato & Hachisu 



[l9m and lKatol (Il991albn . Now a days we have many calculational/observat ional indications 
that fri ctional processe s are ineffective in nova outbursts [ see, e.g., introduction in iKato & Hachisul 
(I2OIII) , Section 7.2 i n Kato. Mikolaiewska & Hac hisu (2012), and also Section \6?2\ of this re- 
view] . Cassisi et al. ( 19981) and jPiersanti et al. I il999> .2000i) extended their arguments on the 
low mass WDs (^ O.BMq) to all the WD masses up to the Chandrasekhar mass, as often cited 
like 'WDs are hardly grow t o reach the Chandrasekhar mass limit due to frictional effects' (e.g. 



Straniero & Piersanti 



20031 lYungelson Il2005l iMaoz & Badenes 11201 01) in the arguments against 



the SD scenario. 



In the SD scenario, on the other hand, we have a much wider mass-growing region of WDs, 
not restricted into the narrow strip in Fig. |2]as in the original DD scenario. In addition to this 
narrow strip, WDs in the region above Mace = M^- can grow in mass during the accretion wind 
evolution. It could also grow in the region slightly below the line M^^c = A^fst, i.e., 10"** M© 
yr"' < Mace < M^t, where the shell flashes are weak. The mass increasing efficiency of WDs is 
obtained from the c ombination (?7H/7He) of the mass accumulation effi ciencies of hydrogen shell 
burning/flashes ( Hachisu et al.lll999bt Prialnik & Kovetz 1995 ) and of helium shell burn- 
ing/flashes T/He (iKato & H ac hisu I I2OO4I) . 



Nelemans I (l2012l) compared the total mass-increasing efficiency of a 1 M© WD between the 
DD and SD models. The mass increasing efficiency is up to rjuriue ~ 0.8 in the SD models, but 
only ?7H?7 He ~ 0. 1 at most ( at log Mace {Mq yr"') ~ -6.6) in the model by Yungelson (2010) (see 
Fig. 4 in iNelemans 2012)). The latter va lue is increased to /yn^/He ~ 0.2 (at log Mace ~ -6.1) 
for symbiotic stars ( Iben & Tutukovlll996 ) because of their larg er Roche lobe sizes. N ow we can 
see where the difference comes from. The original DD model ( Iben & Tutukov ''l996) adopt the 
old opacity and instantaneous envelope mass ejection due to the Roche lobe overflow. Moreover, 
they assumed line-driven winds from the WD of very large mass-loss rate (with 100% efficiency 
of photon momentum flux, i.e., a few times IO^^Mq yr"'). We think that this rate is too large 
and unlikely. Such large mass-loss rates are not included in the SD scenario because no optically 
thick winds blow when Mace < M^r- As a results, the original DD scenario yields a very small 
efficiency of mass-growth of WDs even when they calculate SD binaries in their DD scenarios. 
Thus, the SD channel hardly contributes to the total SN la rate in the DD scenarios. 
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5.6 Other ingredients - accretion disk instability 



King. Rolfe & Schenkea ( 120031) presented an idea that thermal disk instabilities temporarily in- 
crease the mass accretion rate onto the WD, which would trigger steady hydrogen burning, thus, 
the lower line of the narrow strip "steady H-burning" in Fig. |2] would go down as much as by 
-2.4 in logarithmic scale. However, this idea is based on a misunderstanding about the instability 
of hydrogen shel l -burni ng (shell-flash). We should note that the s ame mistake was repeated in 
Alexander et al. ( 201 1 ) after the criticism by iHachisu et al. ( 2010 !). As explained in Section lTTl 
a shell flash occurs when the envelope mass increases to the ignition mass at which the thermal 
balance in the envelope breaks down. This does not occur only when the accretion rate tempor- 
arily exceeds Mg- Unfortunately, this incorr ect idea was alrea dy incorporat ed in some binary 
evolution models dXu & Li 1 120091: IWan feHan 2010a: Meng & Yang 11201 Ol) . but these results 
have no astrophysical meaning. 



6. Related objects 



6.1 Super-Chandrasekhar mass white dwarfs - U Sco, T CrB, and V407 Cyg 



Recently su per-Chandrasekhar mass WDs are suggested from observations of very luminous 
SNe la re.g..lHoweU et al. ''2006: Hi cken et al. ll2007l:[Yamanaka et al. ll2009HScalz o et al. II2OI0I: 
Silverman et al. 201 ll: Taubenberger et al. 201 1 ). In the SD scenario, WDs that e xceed the 



Chandrasekhar mass limit are theoretically possible if they are rapidly rotating (e.g., LFustham 



201 lUDi Stefano et al.l20r It IHachisu et al. l2012aH IHachisu et al .l2012b[ and references therein) 



If recurrent novae are progenitors of SNe la, some of them could be super-Chandrasekhar mass 
WDs . So far, only one ob ject (U Sco) was suggested to harbor a super-Chandrasekhar mass 
WD ( Hachisu et al. Il2012a ). but we here suggest two more objects, the symbiotic recurrent nova 
T CrB and symbiotic classical nova V407 Cyg. 

The WD m ass of U Sco was observa tionally estimated to be Mwd = 1.55 + 0.24 Mq in the 
1999 outburst (Thoroughgood et al. l2001i) . This value is very suggestive, although the ambiguity 
is too large to draw a definite conclusion of super-Chandrasekhar mass. There are additional 
indications of super-Chandrasekhar mass WD in U Sco. (1) The present binary parameters of U 
Sco shows that it just locates in the "final" stage of the path toward super-Chandrasekhar mass 
WDs in the SD scenario (Fig. |3]l, (2) The very fast decline of the optical nova light curve is 
consistent with a very massive WD as large as (or beyond) the Chandrasekhar mass limit of no 
rotation {Mca = 1.4 Mq). (3) The maximum brightness, Mvmax - -7.5, is as much as 3.0 mag 
below the MMRD's prediction (Section [j4b . There is another candidate of super-Chandrasekhar 
mass WD, the symbiotic recurrent nova T CrB. The early decline of its light curve is so rapid and 
its maximum brightness shows a large deviation (2.4 mag) from the MMRD relation (see Section 
13.41 ). We think this kind of large deviations (> 2.5 mag, denoted by asterisk in Tabled from the 
MMRD relation is an indication of super-Chandrasekhar mass WD unless the estimated distance 
is largely different from the true one. 
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V407 Cygni is a long orbital period (D-type) symbiotic star consistin g of a WD and a Mir a 
dMunariet al. Il2012h . V407 Cyg outbursted as a classical nova in 2010 (Ni shivama et al. II2OI0I) . 



The optical light curve shows a very rapid decline in the very early phase which resembles that 
of U Sco. It entered a plateau phase only sev eral da ys after the outburst and started the final 
decline 47 days after the outburst (Munari et al. 1l20l2h . Such an evolution of light curve is faster 
than that of RS Oph. faach isu & KatoJ ( l2012al) modeled the optical light curve of V407 Cyg and 
suggested that the WD mass of V407 Cyg as massive as th at of U Sco and i ts plateau phase can 
be reproduced by a large irradiated disk like in RS Oph (Hac hisu et al. II2OO6I) . The peak absolute 
magnitude is calculated to be Mvniax = -6.6 as in Table |2] and the MMRD error is as large as 
3.0 mag almost similar to that of U Sco. Thus, we expect that V407 Cyg is also a candidate for 
super-Chandrasekhar mass WD. 



Ejecta of the V407 Cyg 2010 nova outburst were enriched with heavy elements dShore et al. 



201 Ih . suggesting that WD core material was dredged up and ejected. This means that the WD 
mass decreases after every nova outburst. Thus, V407 Cyg will possibly not explode as an SN la 
if the WD has a sub-Chandrasekhar mass. If the WD is composed of carbon and oxygen and has 
a super-Chandrasekhar mass, however, it could be a progenitor of SN la because even after many 
nova cycles the W D is still more massive than the Chandrasekhar mass and can explode soon 
after it spins down (IHachisu et al.ll2012bl) . This kind of symbiotic progenitors can explain the 
interact ion between SN la ejecta and circumstellar m atter as seen in Ke pler's supernova re mnant 



mteract ion between aiN la electa ana circumstellar rn atter as seen m Ke pler s supernova re m 
(SNR) jRevnolds et al. Il2007tlchiotellis et al. II2OI2I) and in PTF 1 Ikx dDildav et al.1l2012l) . 



6.2 Helium nova V445 Pup 



V445 Pup is a helium nova that underwent a helium shell flash jKamath & Anupama 2002; 
Kato & Hachisu ll2003l:lAs"hok & Banerjee Il2003l : lwo"udt et al. l[2009l:lKato et aLlbOOS "). The WD 



mass is estimated to be very massive Mwn 



1 .35 Mq, and the companion star is a slightly 



evolved helium star of mass > 0.8 Mp (|Kato et al. I l2008h . The orbital period is suggested to 
be f orb = 0.65 days (Goranskii et al. II2OIOI) . The WD is already as massive as (or close to) the 
Chandrasekhar mass and its mass is incre asing after one cy cle of helium nova. Therefore, this sys- 
tem is a candidate for SN la progenitors (IKato et al ]l2008h . Such a binary, however, does not fit to 
any known path in the binary evolution scenarios (for exampl e, the final or bital periods a re much 



shorter than 0.65 days in their helium star donor channel of IWang et al. 112009: Wang and Han 
201 Obi) , indicating the unknown third path toward SNe la or some missing processes in binary 
evolution models of helium star donor channel. 



In mass increasing WDs like in SSSs or in recurrent novae, we expect that a heUum layer de- 
velops underneath the hydrogen burning zone. This helium layer periodically experiences helium 
shell flashes. V445 Pup indicates that some accreting WDs can become very massive and neither 
helium nor carbon detonation has ever occurred at the bottom of the helium layer. This does 
not support the idea that accreting WDs cannot grow in mass due to frictional process, which is 
frequently assumed in the DD scenarios in order to reject the SD scenarios. 
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